Search for: All records

Single-particle tracking offers detailed information about the motion of molecules in complex environments such as those encountered in live cells, but the interpretation of experimental data is challenging. One of the most powerful tools in the characterization of random processes is the power spectral density. However, because anomalous diffusion processes in complex systems are usually not stationary, the traditional Wiener-Khinchin theorem for the analysis of power spectral densities is invalid. Here, we employ a recently developed tool named aging Wiener-Khinchin theorem to derive the power spectral density of fractional Brownian motion coexisting with a scale-free continuous time random walk, the two most typical anomalous diffusion processes. Using this analysis, we characterize the motion of voltage-gated sodium channels on the surface of hippocampal neurons. Our results show aging where the power spectral density can either increase or decrease with observation time depending on the specific parameters of both underlying processes.

Datasets generated in the report "Aging power spectrum of membrane protein transport and other subordinated random walks". Included data are:

Numerical simulations 
RWdata1.mat: 10,000 realizations, subordinated random walk with Hurst exponent, H=0.3 and \(\alpha\)=0.4.
RWdata3.mat: 10,000 realizations, subordinated random walk with Hurst exponent, H=0.7 and \(\alpha\)=0.4.
RWdata8.mat: 5,000 realizations, subordinated random walk with Hurst exponent, H=0.75 and \(\alpha\)=0.8.
RWdataCTRW.mat: 10,000 realizations, continuous time random walk (CTRW), \(\alpha\)=0.7.

Spectra of simulations
PSDdata1.mat: Power spectral density (PSD) of a subordinated random walk with Hurst exponent, H=0.3 and \(\alpha\)=0.4. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.
PSDdata3.mat: PSD of a subordinated random walk with Hurst exponent, H=0.7 and \(\alpha\)=0.4. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.
PSDdata8.mat: PSD of a subordinated random walk with Hurst exponent, H=0.75 and \(\alpha\)=0.8. Four different realization times are used to compute the PDS: 2^15, 2^16, 2^17, and 2^18.
PSDs_CTRW.mat: PSD of a continuous-time random walk (CTRW), \(\alpha\)=0.7. Five different realization times are used to compute the PDS: 2^8, 2^10, 2^12, 2^14, and 2^16.

Experimental data of Nav1.6 channels in the soma of hippocampal neurons
NavMSDtimes.csv: ensemble-averaged (EA) MSD and time-averaged (TA) MSD. The TA-MSD is measured for three observation times, 64, 128, and 256 frames (3.2, 6.4, and 12.8 s).
NavPSD.csv: Power spectral density (PSD) measured for three observation times, 64, 128, and 256 frames.